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ABSTRACT 

We investigate the 21 cm absorption lines produced by non-linear structures during the early 
stage of reionization, i.e. the starless minihalos and the dwarf galaxies. After a detailed mod- 
elling of their properties, with particular attention to the coupling physics, we determine their 
21 cm absorption line profiles. The infalling gas velocity around minihalos/dwarf galaxies 
strongly affects the line shape, and with the low spin temperatures outside the virial radii of 
the systems, gives rise to horn-like line profiles. The optical depth of a dwarf galaxy is reduced 
for lines of sight penetrating through its HII region, and especially, a large HII region created 
by a dwarf galaxy with higher stellar mass and/or a top-heavy initial mass function results in 
an optical depth trough rather than an absorption line. We compute synthetic spectra of 21 cm 
forest for both high redshift quasars and radio afterglows of gamma ray bursts (GRBs). Even 
with the planned SKA, radio afterglows of most if not all GRBs would still be too dim to 
be the background sources for high resolution (1 kHz) observations, but absorption lines can 
be easily detected towards a high-z quasar. Broadband observation against GRB afterglows 
can also be used to reveal the evolving 21 cm signal from both minihalos and dwarf galaxies 
if there was no X-ray background or it was extremely weak, but it becomes difficult if an 
early X-ray background existed. Hence the 21 cm absorption could be a powerful probe of the 
presence/intensity of the X-ray background and the thermal history of the early universe. 

Key words: line: profiles - galaxies: dwarf - galaxies: high-redshift - cosmology: theory. 



1 INTRODUCTION 

The formation of the earliest galaxies and the cosmic reionization 
are among the milestones in the history of the universe. As the first 
stars form in the earliest non-linear structures, they illuminate the 
ambient intergalactic medium (IGM) and start the reionization pro- 
cess of hydrogen. Based on an instantaneous reionization model, 
the polarization data of cosmic microwave backgroun d (CMB) con- 
strain the redshift of reionization to be z re j on ~ 10.5 l Larson et al.l 
l20ld) . while the Gunn-Peterson troughs dGunn & Petersonll 19651) 
shown in quasar (QSO) absorption spectra suggest that the reion- 
ization of hydr ogen was very nearly complete by z ~ 6 (e.g. 
iFan et al. I l2006h . However, at present we are still unable to ob- 
serve the earliest galaxies directly, and our current understanding 
on the reionization pro cess is only based on theoretical mode ls (e.g. 
iFurlanetto et alj|2004 and simulations (e.g. lTrac et~al]|2008l) . 

The most promising probe of the cosmic reionization is the 
redshifted 21 cm transition of HI which is directly related to the 
neutr al component of the IGM (see e.g. IFurlanetto. Oh & Briggsl 
120061 for a review). Unlike the Lya resonance line, the 21 cm 



E-mail: xuyd@vega.bac.pku.edu.cn 



line could hardly be saturated because of its extremely small Ein- 
stein coefficient (Aio = 2.85 x 10~ 15 s _1 ), so it traces well 
the reionization history, especially during the early stages. Us- 
ing the CMB as the background radio source, 21 cm tomogra- 
phy could map out the three dimensional structure of the emis- 
sion or absorption of 21 cm photons by the IGM (e.g. lMadau et al.l 
1 19971 : iTozzi et alj l200Ch . However, it requires the deviation of 
the spin temperature from the CMB temperature and may not 
be feasible for certain epochs; it could not resolve structures 
smaller than ~ 1 Mpc; and it has a nu mber of observational chal- 
lenges ( IFurlanetto. Oh & Briggsll2006l) . Complementary to the 21 
cm tomography, the 21 cm forest observation detects absorption 
lines of intervening structures against high redshift radio so urces 



jCarilli et alj2002|-|Furlanetto & Loebll2002l : lFurlanetto1l2006t) . and 
it is immune to most of the above difficulties encountered by the 
tomo graphy observa tion. As it is very sensitive to gas tempera- 
ture ixu et alj|2009h . it provides a useful tool to constrain the X- 
ray heating in the early universe. And also, the 21 cm forest is a 
promising probe that could possibly detect high redshift minihalos, 
which are important for determining the mean clumping factor of 
the IGM and putting limits on small scale structure formation. Here 
we focus on the absorption experiments and investigate the 21 cm 
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absorption lines produced by the non-linear structures during the 
epoch of reionization. 

The forest observation relies on the availability of luminous 
radio sources beyond the epoch of reionization. One possibility is 
the high redshift qu asars which have been observed up to z = 6.43 
dWillott eta l. 2007). a nd a r adio-l oud quasar at z = 6.12 was dis- 
covered by McGreer et al.l (120061). Using hi gh-z quasars o r radio 
galaxies as backgrounds. ICarilli et al.l J2002h and lXu et alj J2009h 
have examined the possibility of detecting 21 cm absorption by the 
neutral IGM based on simulations. Another possible option is the 
radio afterglows of high-z gamma ray bursts (GRBs). It is believed 
that some of the long duration bursts are produced by the explo- 
sions of massive stars. The first stars are thought to be likely very 
massive, and may produce bri ght GRBs which would be detectable 
up to a redshift as high as 60 dNaoz & Brom berg 2007 ]). Recently , 
GRB 090423 was disco vered at z = 8.1 dSalvaterra et"ai]|2009l. 
while pTknvir et al .120091 reported its redshift to be z = 8.26), estab- 
lishing the new redshift record of observation for all objects except 
for the CMB. The radio afterglow of a GRB can be observed out 
to a very late time wh en the outflow becomes sub-relativistic (e.g. 
IPihlstrom et al.|[2007h . and considering the cosmic time dilation, it 
offers us an adequate integration time. One possible problem with 
the GRB radio afterglow is that at very low frequency, synchrotron 
self-absorption may become important which may reduce the radio 
flux. Nonetheless, some GRB radio afterglows may be sufficiently 
bright at the relevant frequency ranges. 

In order to plan such observations, we need to know what kind 
of signals will be produced by the early structures, and to under- 
stand the physics behind the expected signals. In this paper, we pro- 
vide a detailed modeling of the 21 cm absorption lines produced by 
minihalos and dwarf galaxies during the epoch of reionization, and 
explore the physical origins of the line profiles. We generate syn- 
thetic spectra of both quasars and GRB afterglows, on top of which 
the 21 cm absorption lines are superposed. Projecting the capability 
of future instruments, we discuss the prospects of detecting the 21 
cm signals from these non-linear objects in the early universe. 

This paper is arranged as follows. In Section [2] we describe 
the physical model involved with the 21 cm absorptions, including 
the halo model for high redshifts, the star burst criterion, a possible 
X-ray background, the physical processes taking place in minihalos 
and dwarf galaxies respectively, and the Lya background produced 
by these early galaxies. In Section[3] we show the spin temperatures 
and 21 cm line profiles of minihalos and dwarf galaxies for various 
parameters. Statistical results including line number density, the- 
oretical spectrum, and the equivalent width (EW) distribution are 
given in Section [4] Then, in Section [5] we study the feasibility of 
the 21 cm observation by making mock spectra of both quasar and 
radio afterglow of GRB. Finally, we summarize and discuss our 
results in Section[6] 

Throughout this paper, we adopt the cosmological parame- 
ters from WMAP5 measurements combined with SN and BAO 
data: Q b = 0.0462, Q c = 0.233, ft A = 0.721, Hp = 
70.1k m s' 1 Mpc"\ <r 8 = 0.817, and n s = 0.96 dKomatsu et all 
120091) . 



2 THE MODEL 

We start with a description of various aspects of physics involved 
with determining the absorption spectrum of 21 cm lines. This in- 
cludes the high redshift number density of dark matter halos, gas 
distribution inside and around halos, the criterion of star forma- 



tion, the X-ray and Lya background, and a detailed modeling of the 
physical properties of starless minihalos and galaxies, respectively. 
Among these properties, the ionization state, the gas temperature 
distribution, and the Lya photon density are especially important 
for determining the strength and line profiles of the 21 cm absorp- 
tion by minihalos/dwarf galaxies. 

2.1 The Halo Model 

In order to model the halo number density at high redshift, we 
use the Sheth-Tormen halo mass function, which is based on an 
ellipsoidal model for perturbation collapse and fits well the simu- 
lation results. The comoving number density of halos at redshift 
z with mass in the in terval (M,M + dM), can be written as 
dSheth&Tormenll 19991) 



n(M, z) dM = F ST (a, z) Q dln ^ . dM 



M dM 



(1) 



where 



F ST (cr,z) = A ■ 



1 + 



ffl<5f c 



exp 



2 a 2 



(2) 

Here po is the cosmic mean density of the total matter today, 
a — cr(M) is the r.m.s. of a Gaussian density field smoothed on a 
mass scale M with a spherical top-hat filter of radius R, where R is 
equivalent to M in a fixed cosmology as R = (3Af/47rpo) 1 ^ 3 , and 
S ac = 1.686/D(z) is the critical overdensity required for spher- 
ical collapse at redshift z, extrapolated to the present time using 
the linear theory, where D(z) is the linear growth factor. The cor- 
rection factors a = 0.707, p = 0.3, and A = 0.3222 were in- 
troduced as appropriate for ellip soidal collapse ( Sheth & Tormenl 
l2002l ; ISheth. Mo & TormerfeoOll) . 

The 21 cm signal depends on the gas content of minihalos or 
galaxies, and the gas fraction in halos will be suppressed by the 
heating processes during the reionization. Therefore, we set the 
lower limit of the halo mass to be the characteristic mass Mc at 
which halos on average could only retain half of their baryons. The 
characteristic mass depends o n the halo merger history and the ther - 
mal evolution of the universe dGnedinll200 0; Okamoto et al ]l2008h . 
At high redshift that we are considering (z > 6), however, the fil- 
tering mass Mf prov ides a good fit to the characteristic mass Mc 
dOkamoto etai1l2008l) . Including the global heating process by an 
X-ray background (see section [231 in the following), we compute 
the thermal evolution of the universe and find that this mass scale 
is ~ 10 6 Mq beyond redshift 7 for an early X-ray background not 
higher than 20% of the intensity today. For high redshifts of inter- 
est, it is not very sensitive to the uncertain intensity of the X-ray 
background because of the delayed response of the gas density dis- 
tribution to the change in the gas temperature. However, if the early 
X-ray background was as high as today, or even higher, it would 
have moderate effect in raising the Jeans mass and the filtering mass 
(see the right panel of FigO. So we take the halo mass range of 
[10 6 M Q , 1O 1O M ] for an X-ray background not higher than 20% 
of today's value, but also consider its effect on the minimum halo 
mass for higher intensities of the X-ray background (see section|4]l. 
The range covers the characteristic halo mass and most of the galax- 
ies th at are responsible for the reionization dChoudhurv & Ferraral 
120071) . 

In the following, we use the NFW density profile for dark 
matter distribution inside of the virial radius r v ir of a halo 
dNavarro. Frenk & Whitelll997h . In general, the key parameter in 
the NFW profile, the concentration parameter c, depends on the 
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halo mass as well as its redshift. Unfortunately, t he concentration 
param eter found from low redshift simulations dCoorav & Shethl 
120020 is not directly applicabl e to the epoch of r eionization. Thanks 
to the resimulation technique, iGao et af] J2005h have simulated an 
example of the first halos from very high redshift with extremely 
high resolution, and its density profile is derived in each resimula- 
tion at the final time. Here we make use of their results. Assuming 
that c is inversely proportional to (1 + z), as indicated by low red- 
shift results, we make a 4th-order polynomial fit to the simulated 
points in the logarithmic space of halo mass. 

For the distribution of gas within the dark matter halo, we as- 
sume that the gas is in hydrostatic equilibrium with the dark matter, 
and have a spherical distribution. We expect gas in the high-z (and 
in general low-mass) galaxies have a rounder distribution rather 
than settling into a disc because: (a) their circular velocities are 
comparable to the sound speed, making the gas more rounded than 
flat; (b) at such high redshifts, there may not be enough time for the 
gas to settle in a smooth disk. With the NFW profile of dark matter 
distribution, one can derive the gas density distribution analytically 
dMakinoetal] 19981) : 



lnp g (r) = lnp g o - 



2 ks T v i r 



[t£(0)-t£(r)], (3) 



where p g o is the central gas density, p is the mean molecular weight 
of the gas depending on the ionization state, m p is the proton mass, 
fee is the Boltzmann constant, and T v i r is the virial temperature of 
the halo. Here v e (r) denotes the gas escape velocity at radius r, 
which can be written as 



v e (r) 



GM(r') 



dr = 2K 2 



F(cx) + 



1 + cx 



xF(c) 



(4) 



where V c 2 = GMj r v i r is the circular velocity at the virial radius, 
x = r/ r'vir, c is the halo concentration, and F(c) — ln(l + c) — 
c/ ( 1 + c) . The maximum escape velocity is reached at the center of 
the halo, u 2 (0) = 2 V c 2 c/F(c). In Eq.©, the central density p g0 
is determined by the condition that the total baryonic mass fraction 
within the virial radius is equal to Sit, /fim, which gives 



Pgo = (Ac/3) c 3 (Q b /Q m ) e A 

Pm f°(l + t) A /H 2 dt ' 



(5) 



where p m = p-m(z) is the mean matter density of the uni- 
verse at redshift z, A = 2c/F(c), and A c is the mean den- 
sity of a virialized halo wit h respect to the cosmic mean value 
d m tervan&Normaij[T998h : A r = 18 tt 2 + 82 (fi^ - 1) - 
39 (p. z m - l) 2 , where Q. z m = Q m (l + zf/ [fi m (l + zf + 0\]. 

To determine the gas distribution outside r v i r , one has to keep 
in mind that the gas slowly falls into the halo because of the gravi- 
tational force of the halo, so the gas has a peculiar (i.e. infalling) ve- 
locity which could be important for the 21 cm line. Also, due to the 
large cross-section at large radii, there is a high probability for the 
line of sight to go through these external parts of t he halos, and this 
infalling gas has to be included in the calculation. iBarkanal d2004h 
has derived a model for the expected profiles of infalling mat- 
ter around viriali zed halos based o n the extended Press-Schechter 
(EPS) formalism teond et aljl99lh . This "Infall Model'Q can also 
be used to compute the final peculiar velocity profiles around viri- 
alized halos as well. The total velocity of the gas is the sum of 
the peculiar velocity and the Hubble expansion. In these infalling 



1 Public code for this "Infall Model" 
http://wise-obs.tau.ac.il/ "barkana/codes.html. 




is available at 



Figure 1. The hydrogen number density profiles in and around dark matter 
halos of different masses: M = 10 6 M Q (solid curves), M = 10 7 M Q 
(dashed curves), and M = 10 s Mq (dotted curves), respectively. The thick 
curves are for redshift 20, while the thin curves are for redshift 10. 



regions around halos, we assume that the gas overdensity and ve- 
locity field follow the dark matter perfectly, with the scaling factor 

f2b/fW 

Several illustrative curves of the total hydrogen number den- 
sity profiles are shown in Fig[T]for various halo masses and two 
sets of redshifts. A constant h ydrogen fraction of primordial value 
Xh = 0.752 is assumed here dSpergel et al .Il2007h . The halos with 
smaller mass and lower redshift are more concentrated. We can see 
that there are discontinuities at r v i r . That is because when we adopt 
the gas density profile above, we are also assuming that the viril- 
ization shock is located right at the virial radius, and at the shock 
radius density and temperature jumps are expected. The exact po- 
sition of the shock radius is subject to debate, and may depend on 
the assembly histor y of the halo. For an cosmology with Q m — 1, 
iBertschingerl dl985l) found that the accretion shock forms at the ra- 
dius R90.8 which enclos es a mean dark ma tter density of 90.8 times 
the cosmic mean, while I Abel et"ai] d2002l) suggests that the shock 
radius is close to the virial radius. Below we assume that the shock 
is located at the virial radius, but note that the exact position of the 
shock radius has only a small effect on the resulting 21 cm signal 
(see the next section). 



2.2 Star Formation Criterion 

After the formation of dark matter halos, some of the halos could 
form stars inside while others could not, depending main ly on the 
cooling processes in each halo (see lMcKee & TarJ 120081 for a re- 
cent review). As for those halos with star formation, a star burst is 
more likely than continuous formation at this early epoch, because, 
small objects are very sensitive to feedbacks and therefore, as soon 
as the first stars formed in the galaxy, their radiation and induced 
supernovae could heat and eject the surrounding gas, quenching 
subsequent star formation. In order to determine whether a halo is 
able to undergo a star burst, here we use a timescale criterion for 
star formation. If the time required for a halo to start forming stars 
is longer than the Hubble time at the halo redshift, it will remain 
starless, i.e. a system that is usually identified as a minihalo. Oth- 
erwise, the halo has enough time to cool and collapse to form stars, 
and becomes a dwarf galaxy. The timescale required for turning on 
a star burst is modeled as the maximum between the free-fall time 
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tfi and the cooling time t coo i, i.e. tss ~ max{ tg, i coo i }, where 
tff = (37r/ 32 Gp) 1//2 . As for the minihalos and dwarf galaxies 
in the early universe, the main coolant is the molecular hydr ogen 
dAbel et alkoollBromm et alj2002l : lBromm & Larsorj2004l) . and 
t coo \ is the H2 cooling time scale. Defining t CO oi = T/T, one can 
obtain ( Teg mark et all 19971) 



^cool 



,200 yr 1 + 



10 T, 



7/2 



60 + T, 4 



512K/T 



(fni) 



(6) 



where T 3 = T/10 3 K, m = n///lcm" 3 , and the / is the fraction 
of hydrogen in molecular form. This fraction varies with time and 
to get the precise value at virialization it would be necessary to 
follow the entire collapse history, but after t he virialization the H 2 
fraction of a halo always saturates at ~ 10 -3 dTegmark et al.ll997n . 
For simplicity, we use / = 10~ 3 hereafter. 

The cooling and collapse processes start from the halo for- 
mation time, which is usually defined as the time at which a halo 
first acquires half of its final mass. There are several predictors 
for the formation redshift di stribution of dark m atter halos, such 
as the standard EPS model jLacev&Colelll993h . the ellipsoidal 
collapse model dSheth & Tormen 2002), and the non-spherical col- 
lapse boundary model 1 Chiueh & Leell20p"l|). At high red shift. the 
EPS model fits the simulation results best ]Lin et alj2003h . Accord- 
ing to the EPS model dLacev&Coiell 19930 , the cumulative proba- 
bility that a halo with mass A/2 at redshift Z2 has a formation time 
fF prior to t\ or a formation redshift zf higher than zi is 



P(t F < ti\M 2 ,t 2 ) = P(z F > z 1 \M 2 ,z 2 ) 



M2 
Mi 



fs 1 {Si,5 sc i\S 2 ,5 3c2 ) dSi, 



(7) 



where fstiSi, 5 3C i\S 2 ,S sc2 ) dSi is the conditional probability 
function: 

Sscl — S sc2 



fs 1 {Si,5 S ci\S 2 ,5 sc2 )dSi = 



2n (Si - S2) 3 / 2 



x exp 



(Sscl — S Bc2 ) 2 



2 (Si - S 2 ) 



dSi (Si > S 2 ,8 ac i > S sc2 ). 



(8) 



Here S is defined as the variance of density field smoothed on a 
mass scale M, Si = o 2 (Mi), S 2 = a 2 (M 2 ), S h = a 2 (M 2 /2), 
S 3C i = Ssc(zi), and S sc2 = S sc (z 2 ). Given any halo mass at the 
redshift of interest, we differentiate Eq.lfTJ with respect to zf to get 
the differential distribution of the halo formation redshift, then its 
formation redshift zf is chosen by Monte-Carlo sampling of the 
distribution curve. 

The star burst time is simply t s = £f + £sb- When the criterion 
t s < tu is satisfied, a star burst occurs at time t s , while t s > tn 
implies that the halo is not able to form stars within a Hubble time. 

2.3 The 21 cm Optical Depth 

With a high redshift GRB afterglow as the background source, the 
neutral hydrogen in the IGM absorbs the 21 cm photons along the 
line of sight. For the diffuse IGM beyond those regions affected 
by the gravity of non-linear structures, it expands uniformly with 
the Hubble flow, then the optical dep th is dFieldll 959: Mada u et al.l 
ll997l:lFurlanetto. Oh & Briggsll2006l) 



TIGM = 



3h P c 3 Aio n m (z) 



32irk B v 2 T s H(z) 



(9) 



where hp is the Planck constant, c is the speed of light, A w = 
2.85 x 10 -15 s _1 is the Einstein coefficient for the spontaneous 



decay of the 21 cm transition, and 1^10 = 1420.4 MHz is the rest 
frame frequency of the 21 cm transition. Here tihi, Ts, and H 
are the neutral hydrogen number density, spin temperature, and the 
Hubble parameter, respectively. 

In order to compute the nm and Ts, a global description of 
the ionization a nd thermal evolution of the IGM is re quired. Sev- 
eral authors (e.g. lOh & Haimar]20 03; Furlanetto 2006) have shown 
that X-rays will likely heat the IGM significantly before reioniza- 
tion completes. Also, if it existed, it could partially ionize the IGM, 
suppress the formation of low mass minihalos, and heat the infall 
regions around surviving minihalos and dwarf galaxies. However, 
the heating processes are quite uncertain during the early stages 
of reionization, and it is unknown what the X-ray intensity should 
be at high redshifts. Here we include only the X-ray background 
heatin g which is the most important heating mechanism dFurlanettol 
l2006h . Calibrating the X-ray (> 0.2 keV) lu minosity of high-z ob- 
jects to nearby starburst galaxies, we have dGilfanov et al]|2004 
lFurlanettoll200d) 



3.4 x 10 40 fx 



SFR 



1 M Q yr- 



ergs 



(10) 



where SFR is the star formation rate, and fx is a correction factor 
accounting for the unknown properties of X-ray emitting sources in 
the early universe. We will set fx = 0.1 as the fiducial value as the 
early galaxies are probably much less productive of X-rays than the 
nearby starburst galaxies, but we take fx as a free parameter and 
vary the level of X-ray background by choosing different values 
of fx to illustrate the sensitivity of our results to this uncertain 
process. Assuming that the star formation rate is proportional t o 
the rate at which matter collapses into galaxies dFurlanettoll200o1) . 
the total X-ray emissivity ex can be written as: 



(x 



_r y1(l %, ( f* d/coll/dz 1 + Z \ 



where n is the total particle number density (including hydrogen 
atoms, protons, electrons, and helium), / co n is the fraction of mat- 
ter in collapsed halos with M > M m i n , where M m i n is typically 
taken to be the mass threshold for atomic cooling which has a virial 
temperature TVir = 10 K, and /* is the star formation efficiency 
normalized to the cosmic baryon content. We use a mass-depe ndent 
/* provided by the handy fit in Salvadori & Ferrari] d2009l) . and 
evaluate it at M m i n for the most populous galaxies for approxima- 
tion. 

A fraction of energy in the X-ray background contributes to 
the heating of the IGM, and another fraction contributes to the ion- 
izatio n of IGM, these fraction s depend on the ionizing state of the 
IGM (Valdes & Ferrara 2008). In the absence of an ionizing UV 
background, the ionization would be caused mainly by these X- 
ray photons, and we take into account both the primary and the 
secondary ionization processes in computing this Xi(z), assuming 
ionization equilibrium. Then th e evolution of the gas temperature 
in the IGM can be described by dFurlanettoll2006h 

^ = -2H(z)T K + 2 -^+ 2 -^, (12) 
at 6 kbti 6 kbti 

where the first item on the right-hand side describes the adiabatic 
cooling due to Hubble expansion, e C om P is the Compton heat- 
ing/cooling rate per unit (physical) volume, and ex,h denotes the 
part of the total X-r ay background g oing into heating. The Comp- 
ton term is (see e.g. lFurlanettoll2006t) 



2 fcomp 

3 fen 



Xi 



H <7T U-/ 

1 + /hc + Xi 3m e c 



(r 7 -T K ; 



(13) 
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Figure 2. The evolution of the global temperature (left panel), the mean ionized fraction {central panel), and the filtering mass and Jeans mass (right panel) of 
the IGM. The dashed curves are for fx = 0, and the solid curves from bottom to top in each panel are for fx = 0.05, fx = 0.1, fx = 0.2, fx = 1.0, and 
fx = 5.0, respectively. In the right panel, the thin and thick curves illustrate the evolutions of the Jeans mass and the filtering mass respectively. 



Here /h c is the He fraction by number, or = 6.65 x 10~ 25 cm 2 is 
the Thomson scattering cross section, u~, = (4a/c)T* oc (1 + z) 4 
is the energy density of the CMB photons, where a = 5.67 x 
10 _5 erg cm -2 s -1 K -4 is the Stefan-Boltzmann constant, and m e 
is the electron mass. The evolutions of the IGM temperature and 
the mean ionized fraction caused by X-rays are shown in the left 
and central panels of Fig[2]for several values of fx- The case of 
fx — is denoted by the dashed curves. It corresponds to the 
situation in which there is no X-ray background, and the IGM tem- 
perature decreases adiabatically with a mean ionized fraction of 
3 x 1(F 4 which is the residual electron fraction left over after the 
recombination. The solid curves from bottom to top in each panel 
take fx = 0.05, fx = 0.1, fx = 0.2, fx = 1, and fx = 5, respec- 
tively. Here we also illustrate the corresponding evolutions of the 
Jeans mass (thin curves) and the filtering mass (thick curves) in the 
right panel of Fig[2] The filtering mass, which provides a reason- 
able fit to the characteristic mass, is M ¥ ~ 10 6 M o for fx < 0.2, 
M F ~2x 10 6 Mq for fx ~ 1, and M F ~ 5 x 10 6 M Q for fx ~ 5 
at z ~ 10. 

The IGM creates a global decrement in the afterglow spec- 
trum, on top of which minihalos and dwarf galaxies produce deep 
and narrow absorption lines. The main broadening mechanism of 
each absorption line is the thermal broadening with the Doppler 
profile. The 21 cm optical depth of an isolated object i s the integral 
of the absorption coefficient along the line of sight dFieldll959l ; 
iMadau et al]l997l ; IFurlanetto & Loebll2002h : 



t(v) = 



3 hp c A A 10 

32 7T 3 / 2 fc B 



nm(r) 
b(r)T s (r] 



exp 



b 2 {r) 



dx, (14) 



where b(r) is the Doppler parameter of the gas, b(r) = 
\J 2 kBTK(r)/mn, u(u) = c(y — fio)/fio, and v(r) is bulk ve- 
locity of gas projected to the line of sight at the radius r. Inside of 
the virial radius, the gas is thermalized, and v(r) — 0, while the gas 
outside the virial radius has a bulk velocity contributed from both 
the infall and the Hubble flow according to the "Infall Model". The 
coordinate x is related to the radius r by r 2 = (a r v i r ) 2 +x 2 , where 



a is the impact parameter of the penetrating line of sight in units of 

The spin temperature of neutral hydrogen is defined by the 
relative occupation numbers of the two hyperfine structure levels, 
and it is determined by three competing processes: (1) absorp- 
tion of CMB photons; (2) collisions with other hydrogen atoms, 
free electrons, and other species; and (3) scattering w ith UV pho- 
tons. The equilibrium spin te mperature is given by dFieid|[l958l ; 



tons. Ih e equilib rium spin te r 
iFurlanetto, Oh & BriggsKOOrjl) 



Ts 1 = 



(15) 



1 ~\~ 3Cc ~\~ 

where T 7 = 2.726(1 + z) K is the CMB temperature at redshift 
z, Tk is the gas kinetic temperature, and Tq is the effective color 
temperature of the UV radiat ion. In most cases, Tq = Tk d ue to 
the frequent Lya scattering jFurlanetto. Oh & Briggj |2006l). The 
collisional coupling is described by the coefficient x c , and x a is 
the coupling coefficient of t he Lya pumping effect known as the 
Wouthuysen-Field coupling jWouthuvser]|l952l; iFieldll 19581) . The 
main contributions to x c are H-H collisions and H-e~ collisions, 
and it can be written as 

„ , OH rp _ HH rp 

eH . HH n e K 10 1 * n H lKio J* ,,,, 

x c = x c + x c = — 1 , (16) 

Aio T 7 Aio 

where = 0.0682 K is the equivalent temperature of the en- 
ergy slitting of the 21 cm transition, and k^q and the 
de-excitation rate coefficients in collisions with free electrons and 
hydrogen atoms respectivel y. These two coefficients at dif ferent 
temperatures are tabulated in IFurlanetto. Oh & Briggsl d2006t) . The 
coupling coefficient x a is proportional to the total scattering rate 
between Lya photons and hydrogen atoms, 

4P„ T, 



27A W T 1 
where the scattering rate P a is given by 



= 4.7T <T a Ja 



(17) 



(18) 



Here a a = f a where f a — 0.4162 is the oscillator strength 

m e c' 
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of the Lya transition, n^ ot is the total number density of Lya pho- 
tons, J a is the number intensity of the Lya photons, and Akd = 
(b/c) v a is the Doppler width with b being the Doppler parameter 
and v a being the Lya frequency. 

In addition to the global Xi{z) and 7igm(z), to compute the 
line profiles of the 21 cm absorptions by minihalos and dwarf 
galaxies, we need a detailed prescription of the ionization state, the 
temperature profile, and the Lya photon density in and around these 
objects. We model respectively these properties of both minihalos 
and dwarf galaxies, as well as the intensity of Lya background in 
the following. 



,x 10 



I 3^ 



2.3.1 Modeling the minihalos 

Minihalos refer to those small halos that are not capable of hosting 
stars, and they represent a very numerous population according to 
the halo mass function. As we are considering the early stages of 
reionization, when the IGM was slightly ionized with some rare 
sites illuminated by the earliest galaxies, the background radiation 
field has not been set up. However, the gas in minihalos, which 
has much higher density, could be collisionally ionized partially, 
depending on its temperature. The HI fraction is computed from 
collisional ionization equilibrium (CIE) 



10 



n e nm 7 = ftB« e n p 



(19) 



where 7 is the collisional ionization coefficient (Cen 1992), and 



ob is the case B recombination coefficient dHui & Gnedinl 



1991) 



which is appropriate for the high-z and low-mass halos as well as 
the IGM that are far from fully ionized. 

The coefficients 7 and as are both functions of temperature. 
During the virialization process, the gas inside a minihalo is shock 
heated to the virial temperature. As for the gas outside the halo 
(r > TVirX since we are looking at high redshifts when reioniza- 
tion is still very patchy, the IGM around a minihalo might be much 
colder or much hotter than the gas inside the minihalo depending 
on its environment. In other words, either the minihalo is embed- 
ded in an almost neutral region, or it is embedded in an HII region 
produced by a neighboring galaxy. If we go to high enough red- 
shifts, the probability for the halo to be in a neutral region increases. 
With an X-ray background, we assume that the temperature of gas 
around each minihalo equals to Tigm • The HI column density of a 
minihalo (integrated from — r V j r to +r v ir along a line of sight pass- 
ing through the halo center) as a function of its mass is shown in 
Fig(3]for redshift 10. For low mass minihalos, the HI column den- 
sity increases with the halo mass, but it decreases rapidly for halo 
masses larger than ~ 5 x 1O 7 M0, as the gas is gradually ionized 
in those halos with higher virial temperatures. 

Another ingredient for the computation of 21 cm optical depth 
is the number density of Lya photons, which is essential for de- 
termining the spin temperature of hydrogen. The background Lya 
photons can not penetrate into minihalos and contributing to the 
Lya flux inside them (see the arguments in section [2.3. 3t . and in 
the absence of stars, the Lya photons could come only from re- 
combinations in the gas itself. The Lya production rate is 



riLya = (2/3) as n e n p . 



(20) 



The factor of two-thirds accounts for the probability of the re- 
combination leading to a Lya photon and the other one-third is 
the prob ability of obtainin g photons of frequencies different from 
the Lya dOsterbrockll 19891) . Because of the extremely large cross- 
section of Lya resonant scattering, the Lya photons cannot escape 
immediately from the halo, but scatter with the hydrogen atoms 



10 
M [M 



7 

solar^ 



10 



Figure 3. The HI column density as a function of minihalo mass at redshift 
10. 



frequently and diffuse in the gas. This process is well described 
as a "random walk". Thus, the Lya photons can accumulate in a 
given region for a diffusion time, idiffu = '?ms/( c 'mfp)> where 
Zrms is the r.m.s. distance that Lya photons have to travel in or- 
der to escape from the region, and Z m f p = (nm a a / Avd)^ 1 is 
the mean free path of the Lya photons. For the case of minihalo 
Irms = TVir- However, the Lya accumulating time f acc is also lim- 
ited by the time interval between the halo formation and the Hubble 
time, i.e. (in — If), then the real accumulating time of Lya pho- 
tons is tacc = min{ tdiffu, 4h — tp }• Actually, tdiffu is quite large, 
and it is always longer than the Hubble time, so the Lya photons 
would accumulate for a time interval (tu — tp). Due to the large 
accumulating time, even with a low ionized fraction in the IGM, the 
Lya photons from recombination could also be important outside 
the r v i r as compared to the background Lya photons. 

In order to illustrate the relative effects of collisional coupling 
and Lya coupling in minihalos, we plot several contours of the co- 
efficients x c (upper panel) and x a (bottom panel) on log(T v i r ) — z 
planes in Fig(4] These coefficients are all calculated for the central 
properties of halos, and to show clearly the dependence of x a on 
the virial temperature of halos, we set the accumulation time of the 
Lya photons to be a Hubble time at each redshift. 

As we see clearly from the contours, x a dominates over x c 
for high mass halos, in which gas is collisionally ionized and Lya 
photons are effectively produced and accumulated from recombi- 
nations. On the other hand, x c dominates over x a for low mass 
halos which are almost completely neutral. 

As for x a , it has a very weak dependence on the redshift, but 
increases dramatically around T v i r ~ 10 4 K due to the change of 
ionized fraction there. At higher temperatures, when the gas is al- 
most collisionally ionized, the x a slowly drops again because of 
the decreasing «b with temperature. 

There are several effects taking place to determine the behav- 
ior of x c . In general, x c increases with temperature because of 
the increasing collisional de-excitation rate k 1q and k^q . At low 
temperatures, x c comes mainly from the H-H collisions, and the 
slowly decreasing x c with increasing redshift is caused by the de- 
creasing central density of hydrogen in less concentrated miniha- 
los at higher redshifts. The coupling changes from H-H collisions- 
dominated to e~-H collisions-dominated as the temperature goes 
up, and this is in part responsible for the non-monotonic behav- 
ior of x c at log(Tvir) ~ 4 — 4.5. In addition, given a redshift, 
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Figure 4. Contours of collisional coupling coefficient x c (upper panel) and 
Lycf coupling coefficient x a (bottom panel) on log(T v j r ) — z planes. From 
bottom to top in each panel, the values of the contours are 7 X 10 3 , 10 4 , 
3 X 10 4 , 5 X 10 4 , 7 X 10 4 , and 10 5 for x c , and 10" 4 , 1, 10 4 , 10 8 , 10 12 , 
10 13 , 10 13 , and 10 12 for x a , respectively. 



the concentration is lower for higher mass (higher T v i r ) halos, and 
the central density of gas is lower. Therefore, n e is smaller at high 
temperatures when the gas is almost ionized and n e represents the 
total hydrogen density. This could also result in the non-monotonic 
behavior of x c . 



2.3.2 Modeling the dwarf galaxies 



As discussed in section [Z2l star bursts could occur in some of the 
halos, turning them into dwarf galaxies when the cooling process 
allows the gas to cool and collapse within a Hubble time. Depend- 
ing on the initial mass function (IMF) and star formation efficiency, 
the stars produce a radiation field that will photonionize the gas cre- 
ating an HII region, heat the gas around, and produce Lya photons 
at the same time. All these effects could have influences on the 
strength and shape of the 21 cm signal of the galaxy. 

The initial mass function and the emission spectrum of the 
first galaxies remain very uncertain (see[ciardi & Ferrara 20051 for 
a review). However, ISchaerea J2002l 120031) has examined spectral 
properties of the ionizing continua of high redshift starburst galax- 
ies for various IMFs, metallicities and star formation historieo We 
make use of their results for the emitting rate of H and He + ionizing 
photons per solar mass of burst stars, denoted by Qh and Q Hc +> as 
well as the average energy per photon with E > 54.4eV (He + 



2 Data are available at: 
http://cdsarc.u-strasbg.fr/cgi-bin/Cat7VI/109 



ionizing threshold) denoted by E Hc + . The Lyman continuum pho- 
tons are divided into two parts: photons with E < lOOeV could 
effectively result in photonionization, and those soft X-rays with 
E > lOOeV have a large probability to escape, and serve as an ex- 
tra heating mechanism in addition to the X-ray background outside 
the HII region. Assuming a simple power law for the high energy 
tail of the emission spectrum, we extrapolate Q Ho + an d to 
get the production rate of soft X-rays Qx and their average pho- 
ton energy Ex- Thus the production rate of ionizing photons is 
Qion = Qh — Qx- 

As we are considering an early reionization epoch when the 
UV background has not been set up, all the ionizing photons cre- 
ating an HII region come from the stellar sources inside of the 
dwarf galaxy, i.e. the Qion. Considering the soft spectra of stel- 
lar sources, we assume a sharp boundary for each HII region, i.e. 
Xi — 1 for r < Run, x% = Xi(z) for r > i?Hii and r > r v i r , and 
Xi — Xi (CIE) for Run < r < rvir when Ran < r v ir. Further, 
we assume all the stars to be located at the center of the galaxy, and 
neglect the Hubble expansion during the growth of the HII region, 
as its growth time is much shorter than the Hubble time. Starting 
from the star burst time, t#, the evolution equation of the radius of 
an HII region is 



d /? 

n m (R)47TR 2 ^ = 7V 7 (i) 



R(t) 



qb n e n p Anr dr, (21) 



where iV 7 (t) is the ionizing rate, and the second term on the right- 
hand side accounts for recombinations in the HII region. Here 
Qb = a B (10 4 K) = 2.59 x 10~ 13 cm 3 s _1 is used, as appro- 
priate for the temperature of an ionized region. Here we negelect 
the change in gas density profile due to the heating, as we will see 
later that it is quite weak. Before the ionizing front breaks out of 
the virial radius, the ionization is caused by ionizing photons that 
did not escape the galaxy, then 



N^t) = (1- /esc) Qion/, 



^ 771 



M. 



(22) 



where / osc is the escape fraction of ionizing photons, and /* is 
the star formation efficiency normalized to the cosmic baryon con- 
tent associated with a halo of mass M. We use f asc = 0.07 
as fav ored by the early reionization model (ERM. iGallerani et al.l 
20081). and a mass-depend ent /* provided by the handy fit in 
Salvadori & Ferrard J2009h . After the ionizing front grows larger 
than the virial radius, all the ionizing photons have escaped, then 
/ csc = 1, and 



JV 7 (t) = /« 



«^ion /* — — M. 



(23) 



We plot the radius of the HII region as a function of the time 
after star burst for a dwarf galaxy with M = 10 7 A/q at redshift 10 
in Fig(5] The (physical) virial radius of this galaxy is r v i r = 641pc. 
Since Qh has a weak dependence on the metallicity but very sen- 
sitive to the IMF, we fix the metallicity at Z = 10~ 7 , and show the 
growth of i?Hii for three di fferent IMFs. These IM F models are 
correspondingly chosen after Schaerer (2002, 2003), with a power- 
law slope of the Salpeter value (a = 2.35) and different upper and 
lower mass limits: 1 - 100 M Q (IMF- A), 1 - 500 M Q (IMF-B), 
and 50 - 500 M e (IMF-C). The Run grows very fast at the begin- 
ning, but limited by the speed of light. At later stages, it is also lim- 
ited by the stellar lifetime. As the stars fade away, the subsequent 
supernova explosions will sweep the surrounding gas away, thus 
suppressing the 21 cm absorption signal. For this reason we stop 
the evolution of the HII radius at its maximum value. As one could 
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Figure 5. The growth of Run for a dwarf galaxy with mass M = 10 7 M Q 
and redshift z = 10 after a star burst. The solid (black), dashed (blue), and 
dotted (red) curves are for IMF-A (1-100 Mq), IMF-B (1 - 500 Mq), and 
IMF-C (50 - 500 Mq), respectively. All the IMF models have a Salpeter 
shape, the metallicity shown here is Z = 10~ 7 , and we set fx = 0.1 for 
the X-ray background. 



expect, the size of an HII region is very sensitive to the IMF due to 
the strong dependence of Lyman continuum production on stellar 
mass. Here we are considering isolated HII regions around dwarf 
galaxies, but in reality the galaxies will be somewhat clustered. So 
the ionized regions may extend further. However, the clustering of 
the first galaxies is beyond the scope of this work, and we reserve 
this to future works. 

Once an HII region is created, the ga s inside rapidl y ap- 
proaches the temperature around 2 x 10 4 K l lMeiksiiJ[2007h . For 
the 21 cm absorption which requires neutral hydrogen, it is more 
important to compute the gas temperature outside the HII region. 
Here we include the heating by local soft X-rays emitted by the 
dwarf galaxies. While the stars in the galaxy are also emitting Lya 
photons, t he heating effect due to the rep eated Lya scattering is 
negligible l ichen & Miralda-Escudell2~004l) . So the evolution of gas 
temperature can still be described by Eq.dl2t, but with an extra con- 
tribution to the ex,h from the local X-rays. The energy deposition 
rate of local soft X-ray heating is 



ex, g = TZnm (Ex — E t h_) /x,h, 



(24) 



where 1Z is the ionizing rate by the soft X-rays, E t h = 13.6eV 
is the ionization threshold of hydrogen, and /x,h is the fraction of 
the primary electron energy deposited into heat. The X-ray ionizing 
rate can be written as TZ = Qx(i) A/* <Ji(v)/(Anr 2 ), where M* = 
/, (Q|,/O m ) M, and ci( v) is the photonionization cross section of 
hydrogen dMeiksinll2007l) : 



Gl(v) 



<?0 



- 



V^th 



V 

'Ah 



(25) 



ing to the ionization threshold 13.6eV, /3 = 1.34, and s = 2.99. 
The average frequency of photons with E ^ 0.1 keV is us ed here. 
As for /x,h, we use a handy fit bv lValdes & Ferraral ( 120081 ), which 
is a function of ionized fraction Xi. 

With the initial condition of the IGM temperature at the star 
burst time, we solve for the gas temperature around an HII region, 
and plot the results for a dwarf galaxy with M = 10 7 Afe at red- 
shift 10 in Fig[6] We take fx = 0.1 as the fiducial value for the 



Figure 6. The temperature profiles of gas around a dwarf galaxy with mass 
M = 10 7 and metallicity Z = 10~ 7 at redshift z = 10. The solid 
(black), dashed (blue), and dotted (red) curves are for Salpeter IMFs with 
mass ranges of 1 - 100 M© (IMF-A), 1 - 500 M© (IMF-B), and 50 - 500 
Mq (IMF-C), respectively. The temperature profile of gas if there is no star 
formation is also plotted with the dot-dashed (green) curve, and a fiducial 
value of fx =0.1 is assumed for the X-ray background. 



X-ray background. The solid, dashed, and dotted curves are for the 
three IMF models used in Fig[5] For illustration, we have assumed 
that the star burst occurred at the same time as the halo formation, 
and the temperatures are all evaluated at t va i uc = tu — tp, where 
the most probable value for zf is used. For comparison, the temper- 
ature profile for a starless minihalo of the same mass and redshift 
is also shown as the dot-dashed line. Inside the HII region, the tem- 
perature is fixed at 2 x 10 4 K, and it drops down very sharply at the 
ionizing front due to the lower heating rate outside, as we could ex- 
pect for stellar sources. The main heating mechanism outside r v i r 
is the X-ray background heating. 

Apart from ionizing photons and soft X-rays, the dwarf galaxy 
is also producing Lya photons which are essential for determining 
the spin temperature of hydrogen. However, those Lya photons di- 
rectly emitted by the galaxy could hardly contribute to the coupling 
in the gas outside the HII region due to the extremely long diffu- 
sion time, and they will be blocked near the ionization front. It is 
the Lya photons cascaded from soft X-rays that penetrate into the 
nearby IGM and effectively couple the spin temperature to the ki- 
netic temperature of the gas l lChen & Miralda-Escudel2008l) . So the 
contributing Lya photons near an HII region come from recombi- 
nation, the Lya backgro und, as well as the soft X-ray cascading. 
IValdes & Ferraral d2008l) have examined the fraction fhya of the 
primary energy of electrons ionized by X-rays that converts into 
Lya radiation. Using their result, the Lya production rate from cas- 
cading around the HII region is 



t-, Ex — E t h , 

"■Lya = K-rim - /Lya, 

hV a 



(26) 



where v a = 2.47 x 10 15 Hz is the frequency of Lya photons. 
The total number density of Lya photons is obtained by integrating 
riLya over an accumulating time i acc . 

In this case, the accumulating time of Lya photons is t &cc = 
min{ tdiffu, £xray, (ts—ts) }, where tdiffu is the Lya diffusion time 
as before, but we take the mean free path (m.f.p.) of the soft X-rays 
as the i rms here, fxray is the time during which the dwarf galaxy 
is producing X-rays, and (tu — t s ) is the time interval between the 
halo redshift and the star burst redshift. Because of the large m.f.p. 
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of the X-rays (compared to the r v ir), the Lya diffusion time is also 
very large in this case, even larger than the Hubble time. As a result, 
the Lya accumulating time is determined by (tu — t s ) for the cases 
with IMF-A and IMF-B, while for the case with IMF-C, the Lya 
accumulating time is limited by the lifetime of the stars. 

2.3.3 The Lya background 

In addition to the Lya photons produced by recombination and 
those cascaded Lya photons near galaxies, there is also a Lya 
background. The Lya background originates from the continuum 
photons to the blue side of the Lya line emitted by the early 
galaxies, and its intensity beco mes important at very high redshift 
Chen & Miralda-Escu de 2008 ). Only photons with v < vp can be 
redshifted to v a , and all photons of higher frequencies will be ab- 
sorbed at the Ly/3 or higher resonance Lyman series lines. Let E(y) 
be the number of photons emitted over stars' lifetime per baryon in 
one of the halos that collapse to form stars and per u nit of fre- 
quenc y. According to the stellar spectra we adopted from lSchaerej 
J2003h . E(v) is well approximated as flat near the Lya frequency, 
the n the intensity of the continuum Lya background can be written 
as l lChen & Miralda-Escudell2008l) 

Ja = Jo E a v a [F{z) - F{z p )]. (27) 

Here Jo is the fiducial Lya intensity corresponding to a photon 
density of one photon per hydrogen atom per log frequency, 



where tih is the number density of hydrogen. E a — E(v a ) is 
averaged over all the galaxies in the mass range of our considera- 
tio n, and we estim ate this value from the spectrum data provided 
bv lSchaerej d2003h . In Eq, J27b . F(z) is the fraction of mass bound 
in star-forming halos, 

fM n , ax M 

F(z) = / — n(M,z)p i ,(M,z)dM, (29) 
JM min Po 

where n(M,z) is the halo mass function and p*(M,z) is the 
probability of having star formation for a halo with mass M 
at redshift z. We derive p*(M, z) according to the star forma- 
tion criterion given in Section 12.21 At z = 20, J a ~ 3 x 
10" 11 cm" 2 s^sr" 1 Hz -1 , while at z = 10, it has increased to 
J a ~ 2 x 10~ 10 cm" 2 s _1 sr _1 Hz -1 which is already important 
for the Wouthuysen-Field coupling effect. 

The problem of Lya background propagating into minihalos 
is complicated. The situation in minihalos would be different from 
that in the IGM because of their higher column densities of neutral 
hydrogen (typically, Nm ~ 10 21 cm" 2 ). The Lya scattering cross 
section in a region of high density consists of a Doppler core with 
a typical width of A1/0 = b/X a centered on the Lya frequency, 
and Lorentzian wings outside the core. Once a photon is redshifted 
or scattered into the "core" frequencies, the cross section is large 
and its mean free path becomes very short and is spatially confined 
within a small region with very long diffusion time. Thus, the pho- 
tons which eventually make up the Lya background come from the 
blue side of the Lya line, at the edge of the "core" and "wing". 
In the case of minihalo, however, these photons would be blocked 
at the surface of the minihalo. The redshift across the minihalo is 
also too small for bluer wing photons which is not blocked at the 
edge of the minihalo to get redshifted into the core-wing boundary 
region when they arrive at the center of the halo. 

Another possibility is to consider photons with such frequency 



that their optical depth through the minihalo is about 1 . Such a pho- 
ton would scatter once inside the halo. One scattering per photon 
is of course negligible for the coupling effect in itself, but the fre- 
quency of the photon would be changed after the scattering. If the 
frequency increases, then the photon cross section become even 
smaller and would escape right away, but if the photon loses en- 
ergy during the collision and shifts to lower frequency, the cross 
section would be increased. Could such a change bring the photon 
from the "wing" to the "core"? The frequency v of such a wing 
photon is given by Nm{r v i r ) cr a (v) = 1, where Nm{r vir) is the 
column density of the minihalo from r v i r to its center, and a a (y) is 
the scattering cross section on the Lya damping wing. For typical 
thermal velocity of the H atom in the minihalo, we estimate that the 
change of photon frequency in one scattering d v < 2.2 x 10 Hz, 
comparable to the size of the core. On the other hand, the frequency 
distance to the line center for such a photon is {y — v a ) ~ 10 13 Hz, 
so df ~ \%(y — Va). Therefore, during one or even a few fa- 
vorable scatterings the photon could not enter the core and should 
escape away. We conclude that with the exception of the surface, 
the Lya background could not affect the spin temperatures of the 
minihalos. 



3 THE 21 CM LINE PROFILES 

With the detailed models of gas density, velocity profile, ioniza- 
tion state, temperature evolution, and the intensity of both local and 
global Lya photons for the minihalos and the dwarf galaxies de- 
scribed above, we are now ready to compute their spin temperature 
profiles and 21 cm optical depths. In this section, we show our re- 
sults of these detailed profiles for a variety of parameters. Although 
it will be very challenging to resolve them with radio instruments in 
the near future, these analyses help us explore the physical origins 
of the line profiles and understand the physics behind. 

3.1 The Coupling Effects and Spin Temperature 

In Fig[7] we plot the spin temperature as a function of distance from 
the halo center for a minihalo (left panel) and a dwarf galaxy (right 
panel) of the same mass M — 10 7 A/© at redshift z = 10. We have 
taken fx — 0.1 as the fiducial value for the X-ray background. In 
order to show the coupling effects, the kinetic temperature of the 
gas and the CMB temperature are also plotted. As we see clearly 
from the left panel, the spin temperature of the minihalo is coupled 
to the gas kinetic temperature (T v j r of the halo) at the center. The 
halo with M = 1O 7 M and z = 10 has T vir ~ 5000 K, and 
the gas is almost neutral in collisional ionization equilibrium. As 
a result, the Lya photons from recombination is totally negligible 
and the collisional coupling dominates. As the radius increases, the 
collisional coupling becomes less and less effective because of the 
decreasing density, and the spin temperature gradually decouples 
from the kinetic temperature. In addition to the dropping density, 
the collisional de-excitation rate coefficients (k^} and k±o) also 
decrease due to the sharply decreasing kinetic temperature outside 
r v ; r , dramatically reducing the collisional coupling effect. How- 
ever, the Lya background at z = 10 is already strong enough to 
couple the spin temperature closely to the kinetic temperature of 
the gas. And the recombination in the IGM also produces a signif- 
icant amount of Lya photons. As a result, T3 is always coupled to 
Tk outside minihalos at 2 = 10. 

As for the spin temperature profile of the dwarf galaxy in 
the right panel of Fig[7] the IMF model-A and a metallicity of 
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Figure 7. Spin temperature (dashed lines) and kinetic temperature (dot-dashed lines) profiles of a minihalo/dwarf galaxy with mass M = 1O 7 M atredshift 
10. A fiducial value of /x = 0.1 is assumed for the X-ray background. The solid curves represent the CMB temperature. Left: the case for a minihalo in CIE. 
Right: the case for a dwarf galaxy which is photonionized after a star burst with IMF model-A and metallicity Z = 10~ 7 . The temperature curves are cut at 
the HII radius of 7.02 r v i r . 




r / r . 
vi r 



Figure 8. Spin temperature profiles of minihalos/dwarf galaxies of different masses: M = 10 6 M© (solid curves), M = 10 7 Mq (dashed curves), and 
M = 10 8 Mq (dot-dashed curves), respectively. A fiducial value of fx = 0.1 is assumed for the X-ray background. The set of thick curves is for redshift 20, 
and the set of thin curves is for redshift 10. Left: the case for minihalos in CIE. Right: the case for dwarf galaxies that are photonionized after a star burst with 
IMF model-A and metallicity Z = 10~ 7 . 



Z — 10~ are assumed. We stop the Run and Tk evolution at 
ivaiuc = tu — tp, where 4f is the most probable formation time 
of this halo, and the spin temperature is also evaluated at this time. 
The left-hand side cut-offs of these curves are just the position of 
-Rhii (7.02 r v ir). As the HII radius of this galaxy is larger than the 
virial radius, the low density and low temperature outside the HII 
region make the collisional coupling very weak. In this case, the 
Lya pumping is always the dominating coupling effect. The global 
Lya photons from recombination and background flux dominate 
over the local Lya photons from soft X-ray cascading, and just as 
the case outside a minihalo, the spin temperature sticks to the gas 
kinetic temperature. 

Then we investigate how the spin temperature changes with 
halo mass and redshift, and plot two sets of curves for various halo 
masses in Fig(8] one (the thin set) for z = 10, and the other (the 
thick set) for z — 20. A fiducial value of fx = 0.1 is assumed here. 
For the minihalos in the left panel, the spin temperature is closely 
coupled to the kinetic temperature outside the r v ir for halos at z = 



10, in part because of the Lya background, and in part because of 
the accumulated Lya photons from recombination in the IGM that 
is partially ionized by the X-ray background. But at z — 20, both 
the Lya background and the X-ray background are still weak, and 
Ts lies between the Tomb and Tk for gas around minihalos at this 
redshift. Inside the halos, the density is larger at higher redshift, and 
the collisional coupling effect is correspondingly stronger, so the 
spin temperature is more tightly coupled to the virial temperature. 
On the other hand, halos are more concentrated at lower redshift, 
and the slope of gas density profile is steeper. As a result the spin 
temperature drops more rapidly with radius at lower redshift. One 
different feature for the 10 Mg halo is that its spin temperature 
is effectively coupled to the gas kinetic temperature out to r v i r , 
because for this relatively high mass halo, its virial temperature 
(~ 1.5 x 10 4 K for z = 10, and ~ 2.3 x 10 4 K for z = 20) is 
high enough to make the gas collisionally ionized (xi ~ 52% for 
z = 10, and Xi ~ 98% for z = 20). With this partially ionized 
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Figure 9. Optical depth profiles of a minihalo/dwarf galaxy with mass M = 10 7 M e at redshift 10. All the curves take fx = 0.1. The impact parameters 
of the thick curves are a = (solid cyan), 0.3 (short-dashed blue), 1 (dotted red), 3 (dot-dashed green), and 10 (long-dashed magenta), respectively, and the 
thin solid black curve is for a = 30. Left: the case for a minihalo in CIE. Right: the case for a dwarf galaxy which is photonionized after a star burst with IMF 
model- A and metallicity Z = 10 -7 . 



gas, Lya photons from recombinations are effectively trapped, and 
serve as a strong coupling agent. 

In the right panel of Fig[8] the spin temperature profiles for 
dwarf galaxies are all cut off at their HII radii at the left-hand side. 
For the case of 10 6 Mq dwarf galaxy, due to the much lower star 
formation efficiency (/* ~ 4x 10~ 5 , which increase s to ~ 0.03 for 
a galaxy with 1O 8 M ) JSalvadori & Ferrardl20o"9|) , it could only 
create a very small HII region (i?Hii ~ 0.008 r vlI for z = 10, and 
Run ~ 0.02 r vi r for z = 20). As for the galaxies with 10 7 M Q 
and 10 s Mq, they have larger stellar masses and create larger HII 
regions. So they have much less gas which could contribute to the 
21 cm line absorption. 

3.2 The Absorption Line Profiles 

Consider an isolated minihalo with mass M = 10 7 M© at redshift 
z = 10 for which we compute the optical depth as a function of fre- 
quency using Eq.d 14t. Results are shown in the left panel of Fig|9] 
for different impact parameters of the lines of sight. Very interest- 
ing features can be seen in the shown profiles. First, there are two 
peaks which sandwiched the center of the line for impact parame- 
ters a = 1 and a = 3, showing horn-like profiles; second, as the 
impact parameter increases, the peak optical depth is NOT always 
decreasing. These two points are related issues, and they are both 
caused by the infalling gas with low spin temperature outside the 
minihalo. 

To see clearly the origin of these interesting profiles, we have 
to understand the different contributions to the 21 cm absorption 
from the gas located at different radii. Because of the infall, the 
absorption line produced by the gas at the far side of the halo is 
blue-shifted, while that produced by the gas at the near side of the 
halo is redshifted. Let u p (r) be the peak frequency of optical depth 
created by gas located at radius r, and according to the integrand in 
Eq.d 14t. the absorption is shifted from the line center by 

v(t} 

v p (r) - no = uio. (30) 

c 

Considering a line of sight passing through a minihalo with M = 
10 7 Mq and z = 10 from the center (a = 0), we divide the integra- 
tion into segments, each of which has a contribution from a length 



of 1 TVir on the line of sight. Then we plot every absorption line 
created by one segment in Fig[l0] In the upper panel, the absorp- 
tion line produced by the gas inside of r v j r is shown as the dashed 
curve, the solid lines from right to left correspond to the absorptions 
by segments of (1 — 2) r v j r , (2 — 3) r v i r , (3 — 4) r v i r , (4 — 5) r v i r , 
(5 — 6) r v ir, and (6 — 7) r v i r , respectively, and the dotted-dashed 
line represents the absorption by the segment of (7 — 8) r v tr. In the 
bottom panel, the 12 solid lines from left to right correspond to the 
absorptions by segments of 1 r v i r each starting from 8 r v - a , and the 
last curve on the right represents the integral absorption from 20 to 
100r vir . 

As seen clearly from Fig[l0] the virialized gas inside the mini- 
halo has no bulk velocity, and the peak optical depth is located at 
the center. The gas at (1 — 2) r v i r has the highest infall velocity, 
and the corresponding profile lies at the largest distance to the line 
center in the upper panel. As the radius increases, the t„ profile 
gets closer to the line center because of the lower infall velocity, 
and the optical depth decreases slowly with the decreasing density. 
One interesting feature is that there is a special position where the 
total velocity of the gas changes from negative (infall dominated) 
to positive (Hubble flow dominated), and at this turning point the 
two absorption lines created by the two segments on both sides of 
the minihalo converge into one at the line center, and they con- 
tribute substantially to the central optical depth. For the minihalo 
with M — 10 7 M© and z = 10, this turning point lies at 7.3 r v i r , 
and consequently the absorption by the segment of (7 — 8) r v i r 
peaks at the line center. After that, the r„ profiles that come from 
larger radii leave the line center again, because the infall velocity 
becomes even smaller and the total velocity (in the same direction 
as Hubble flow) becomes more and more positive. Finally, when 
it goes out of the region influenced by the minihalo's gravity, and 
the density drops to the cosmic mean value, we recover the IGM 
optical depth. 

The line profiles in Fig(9]are better understood by noting that 
each is an integral of contributions from different radii. A substan- 
tial contribution to the optical depth comes from the outer region, 
because the gas outside r v ir has lower spin temperature. Especially, 
the gas in the infalling region on the far (near) side of the halo 
shares the same bulk velocity with the gas out of the region on 
the near (far) side, and they absorb the 21 cm photons at the same 
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Figure 10. Different contributions to the optical depth from the gas at dif- 
ferent radii. Upper panel: the absorption by gas inside of r v - lT is shown as 
the dashed line, the solid lines from right to left correspond to the absorp- 
tions by segments of (1 — 2) r v i r , (2 — 3) r v j r , (3 — 4) r vlr , (4 — 5) r vlT , 
(5 — 6) r v ; r , and (6 — 7) r v ; r , respectively, and the dotted-dashed line rep- 
resents the absorption by the segment of (7 — 8) r v - lr . Bottom panel: the 
12 solid lines from left to right correspond to the absorptions by segments 
of 1 r v ; r each starting from 8 r v ; r , and the last curve represents the integral 
absorption from 20 to 100 r v i r . 



frequency. Therefore, the optical depth at the frequency range cor- 
responding to the infalling region is increased a lot. In addition, the 
line profile gets narrower with lower temperature compared to the 
gas inside the halo, which further increases the peak value. As the 
impact parameter a increases from to 1 , more contribution comes 
from this cold region, so the peak optical depth increases. Also, the 
infall velocity shifts the peaks away from the center, and results in 
the horn-like profile. As a increases further up to 3, the infall ve- 
locity decreases, and the two peaks move closer to each other. For 
a larger than the radius of the velocity turning point, the two peaks 
merge together and decrease slowly to the IGM optical depth as 
a — > co. 

The 21 cm profiles for an isolated dwarf galaxy for various 
impact parameters are shown in the right panel of Fig[9] Just as 
before, the spin temperature and the optical depth profiles are all 
evaluated at the time £ va i ue = tu — tp, and the most probable 
value for zf is used. Thus we get an upper limit on the ioniza- 
tion, heating, and Lya coupling effects. The line profiles are com- 
pletely different from the case of minihalos. The horn-like profiles 
disappears and the absorption is strongly reduced for small impact 
parameters. The dwarf galaxy with M = 1O 7 M0 and z = 10 
has an HII radius of 7.3 r v j r , which is close to the turning point 
of the gas velocity. In other words, the hydrogen atoms inside r v j r 
and those within the infalling region are totally ionized, erasing the 
absorption features contributed by hydrogen in this region. There- 
fore, there will be only one peak at the center no matter what the 
impact parameter is, and the optical depth is reduced for lines of 
sight which penetrate the HII region. As the impact parameter a 
increases, the optical depth first increases because more neutral gas 
near the HII region is intercepted by the line of sight. It reaches a 



maximum when a ~ i?Hii and then drops, approaching the IGM 
optical depth. 

We also show the optical depth profiles for different halo 
masses and redshifts in FigQT] Here we also take fx — 0.1 for 
the X-ray background. In the left panel, the case for minihalos is 
shown. As we expected, the line profiles are broader for halos with 
greater masses. That is because, on one hand, the virial temperature 
is higher for halos of greater masses, which results in a broader 
Doppler profile for the absorption by gas inside the halo, and on 
the other hand, more massive halos have stronger gravitational in- 
fluence on the surrounding gas, and the induced higher infalling 
velocity shifts the absorption line farther from the line center. In 
general, the absorption is stronger at higher redshift. The reason is 
that the IGM is denser at higher redshift, and also, the X-ray back- 
ground is gradually set up as redshift decreases, and it heats and 
partially ionizes the gas in the IGM. As the frequency gets far from 
the line center, the optical depths of all the minihalos approach to 
the mean IGM value. 

In the right panel of Fig[TT] we show the 21 cm line profiles 
for dwarf galaxies with the same masses and redshifts as those 
for minihalos in the left panel. Star bursts with IMF model-A and 
metallicity Z — 10~ 7 are assumed. An interesting feature emerges 
for the galaxy with 10 s Mq. Because of the higher star forma- 
tion efficiency associated with higher mass galaxies, this dwarf 
galaxy creates a large HII region (Run ~ 38r v i r for z = 10, 
and Rmi ~ 34 r v j r for z = 20) erasing all the absorption inside of 
it. As a result, a broad optical depth trough is produced instead of 
an absorption line! In the case of IMF model-C, the dwarf galaxy 
of the same mass could ionize an even larger HII region and hence 
could result in an even broader optical depth trough. 

All the line profiles above are computed assuming that the 
minihalo or the dwarf galaxy is isolated. In real cosmic structures, 
a halo is surrounded by other halos, and if we integrate the optical 
depth to a distance larger than the mean separation D of the halos, 
we will probably hit another halo. So the integration should only 
be considered as reliable up to a distance of D/2. In practice, a 
10 s Mq halo, for example, might have many smaller halos closer 
to it than another 10 8 Mq halo. Therefore, it is the mean separa- 
tion of the smallest halos we are considering, i.e. the halos with 
M = 1O 6 M (for fx < 0.2), that determines the integration limit. 
We denote this mean separation as D m i n . The maximum impact pa- 
rameter of a line of sight should also be a max = Dmm/2. Integrat- 
ing the optical depth up to D m i n /2, we plot optical depth profiles 
of a minihalo (blue dashed curve) and a dwarf galaxy (black solid 
curve) with M = 1O 7 M at z = 10 in Fig[l2] The line of sight 
is assumed to be passing through the minihalo/dwarf galaxy from 
the center, and we set fx — 0.1. We see that the dwarf galaxy only 
produces a narrow and weak absorption line at the center, because 
for this galaxy, only the gas in a sphere between the HII radius 
(~ 7.3 r v i r ) and D m i n /2 (~ 9.1 r vlI ) contributes to the absorption. 
In reality, however, the optical depth will not drop to zero at the 
boundaries of the absorption line but will connect with another line 
created by a neighboring halo. In addition, as we mentioned be- 
fore, the clustering of dwarf galaxies will extend their HII regions. 
Similarly, some minihalos will be clustered around the dwarfs, and 
the surrounding gas in the infall region may be ionized, even if 
the minihalos themselves can self-shield. Therefore, the clustering 
can reduce the 21 cm optical depth of some minihalos and dwarf 
galaxies. Although the problem should ideally include the cluster- 
ing properties of early galaxies, it is beyond the scope of this paper 
to include such features. 
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Figure 11. Optical depth profiles of minihalos/dwarf galaxies with different masses: M = 10 6 M@ (solid curves), M = 10 7 A/q (dashed curves), and 
M = 10 8 Mq (dotted curves), respectively. The thick curves are for redshift 20, while the thin curves are for redshift 10. All the curves take fx = 0.1. Left: 
the case for minihalos in CIE. Right: the case for dwarf galaxies that are photonionized after a star burst with IMF model-A and metallicity Z = 10~ 7 . The 
impact parameter shown here is a = 0. 
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Figure 12. The optical depth profiles of a minihalo (blue dashed curve) 
and a dwarf galaxy (black solid curve) with M = 1O 7 M0 at z = 10. A 
star burst with IMF model-A and metallicity Z = 10 -7 is assumed for 
this galaxy. The optical depth at each point is integrated up to D mln /2, the 
impact parameter is a = 0, and the X-ray parameter is fx = 0.1. 



Figure 13. Theoretical evolution of line number density produced by mini- 
halos/dwarf galaxies along a line of sight. The maximum halo mass is 
1O 1O M , and the minimum halo masses are 1O 6 M (solid), 2 X 10 6 M o 
(dashed) and 5 X 10 6 M Q (dotted), respectively. 



4 THEORETICAL SPECTRUM 

In order to superimpose the 21 cm absorption lines of minihalos and 
dwarf galaxies onto a radio spectrum of GRB afterglow, we first 
have to compute the line number density per unit redshift. Theoret- 
ically, including all (weak and strong) absorption lines, the number 
of halo intersections along a line of sight per redshift interval is 

^ = (1 + Z ) 2 J / ,AW n(M,z)A m ^dM, (31) 

dz dz JM min 

where dr / dz is the comoving radial distance per redshift interval, 
n(M, z) is the halo mass function given by Eq.(|TJ, and A max = 
ira^ax is the cross-section (in physical coordinates) of a halo, in 
which the maximum impact parameter a max is set by half of the 
mean halo separation, i.e. _D m i n /2, at each redshift. We plot the 
line number density as a function of redshift in Fig[T3]for three 
different values of M m - m , which are appropriate for different levels 
of the X-ray background. The curves are cut off at redshift z — 8 



since our model applies only to the early stages of reionization, 
when the stellar sources have not set up an ionizing background. 

This line number density is also the probability P(z) that a 
line of sight intersects an object in a r edshift interval dz centere d 
on z. In analogy to the method used in iFurlanetto & Loebl J2002h . 
we divide the observed frequency band into small bins, each of 
which corresponds to a redshift interval Az that is small enough to 
make sure that the probability to have an intersection in each bin 
P(z) Az < 0.01. Then we generate a random number Ri (uni- 
formly distributed in [0, 1)) for each bin, and an intersection is said 
to take place if the condition Ri < P(z) Az is satisfied. 

When an intersection occurs, we randomly choose a mass M 
for the intersected halo according to the mass function. With the 
halo mass and the redshift of that bin, we compute its formation 
redshift distribution, and zf is Monte-Carlo sampled from the dis- 
tribution curve. Then we use the star formation criterion described 
in section |2~2l to determine whether it is a minihalo in collisional 
ionization equilibrium or a dwarf galaxy photonionized after a star 
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burst. For every dwarf galaxy, it s metallicity is fitted as a fun ction of 
halo mass to the results given bv lSalvadori & Ferra ra ( 2009 ). As the 
fitted metallicity is always > 10~ , the IMF model-A is used. With 
the impact parameter randomly selected with an equal probability 
per unit area for each object, we compute the 21 cm line profile for 
every minihalo and dwarf galaxy intersected by a sightline. From 
129 MHz (z ~ 10) to 158 MHz (z ~ 8), we found 3241 lines if we 
set M m i n = 10 6 Mq, out of which 210 lines are attributed to dwarf 
galaxies. That is about 6.5% of the lines coming from dwarf galax- 
ies, and the rest come from minihalos. The relative transmission 
T = exp(— r) for a range of observed frequency corresponding 
to z ~ 10 is shown in Fig[l4] The six panels from top to bottom 
show the results for fx — 0, 0.05, 0.1, 0.2, 1 and 5, respectively. 
M m in = 10 6 Af o is used for fx < 0.2, M min = 2 x 1O 6 M is 
used for fx = 1, and we use Mmin = 5 x 10 6 Mq for the case of 
fx = 5. Note that the y-axes are different between panels. 

As seen from the figure, the absorption lines are very narrow 
and closely spaced, resembling a 21 cm forest. The depth of the 
21 cm absorption strongly depends on the existence and the inten- 
sity of the uncertain X-ray background. The absorption strength is 
significantly increased if there is no X-ray background, and it de- 
creases rapidly with the increase of the X-ray background, which 
ionizes and heats the gas in the IGM reducing the optical depth of 
both non-linear structures and the global absorption. Especially, if 
the X-ray background in the early universe was very strong (e.g. 
fx = 5), resulted from a large number of high mass X-ray binaries 
in the case of top-heavy IMF towards high- 2 or from mini-quasars, 
we could hardly see any feature on a spectrum. This is because, 
on one hand, the line number density is reduced, and on the other 
hand, the absorption by non-linear structures is so weak that most 
lines are covered by the global absorption. For some minihalos, 
the gas outside r v i r contributes to the absorption far from the line 
center because of the infall velocity, and a few absorption lines of 
the neighboring halos overlap with each other on the wings. In this 
case the actual optical depth of each point is the sum over all the 
absorptions caused by these overlapping lines. 

We compared two spectra: one taking into account the local 
X-rays contributing to heating and cascading to Lya photons, the 
second neglecting the above. It turns out that for the IMF model-A 
and the relatively high metallicities that are used in the spectrum, 
the small amount of soft X-rays produced by the dwarf galaxies 
make little difference (less than 0. 1%) even without an X-ray back- 
ground. So we could not see any signature that can be attributed to 
them in the spectrum. However, if a dwarf galaxy with relatively 
high mass and (almost) primordial composition d oes exist, a top- 
heavy IMF (like the IMF model-C) is expected dSchneider et al.l 
12003, l2003h . and the associated large HII region will result in a 
broad bump on the spectrum, erasing several absorption lines pro- 
duced by neighboring minihalos. 

Directly from the spectrum, we could compute the distribu- 
tion of equivalent width of the absorption lines for a specific range 
of observed frequency corresponding to a specific redshift. As the 
continuum of a background source has a global decrement due to 
the absorption of the diffuse IGM, the real signal of non-linear 
structures is the extra absorption with respect to the flux transmitted 
through the IGM. Therefore, the equivalent width of an absorption 
line should be defined as 
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Figure 14. Relative transmission along a line of sight at redshift z ~ 10. 
The six panels from top to bottom show the results for fx = 0, fx = 0.05, 
/ x = o.l, fx = 0.2, fx = 1.0, and fx = 5.0, respectively. Note that the 
y-axes are different between panels. 



where f c is the continuum flux of the background radio source, 
and tigm(z) is the optical depth of the diffuse IGM at redshift z. 
Using the theoretical spectrum of 129 - 133 MHz (corresponds to 
z = 10.01 — 9.68), we compute the cumulative distribution of 
equivalent width of those 21 cm lines around 2 ~ 10, which is 
shown in Fig[13] The dot-dashed, short dashed, solid, dotted, and 
long dashed curves take fx = 0, 0.05, 0.1, 0.2, and 1, respectively. 
As for the case of fx — 5, the number of lines emerged from 
the IGM absorption in this frequency range is too small to derive 
any statistical implication, so the distribution of equivalent width 
is not shown for this case. For the fiducial value fx = 0.1, the 
majority of absorption lines have equivalent widths around 0.03 to 
0.3 kHz. We see that the number of large signals is very sensitive to 
the presence and intensity of the X-ray background. If there is no 
X-ray background (fx — 0), then about 95% of absorption lines 
have W v > 0.1 kHz. For the value fx — 0.05, there are about 
75% of lines with W v > 0.1 kHz, while for fx = 0.1 and 0.2, 
this fraction drops to 45% and 15% respectively. If fx = 1 or even 
higher, there will be no absorption line gets W v > 0.1 kHz. 



5 THE OBSERVABILITY 

There are two kinds of radio sources which are potentially usable 
for absorption line studies at high redshifts during the epoch of 
reionization. One is the high redshift quasars, which are quite lu- 
minous, but so far no quasar has been confirmed at z > 6.5. The 
other candidate is the radio afterglows of GRBs. Although they are 
fainter than quasars, they are more likely to exist at higher redshifts. 
Also, they have simpler power-law spectra at low frequencies due 
to synchrotron self-absorption, so it may be easier to extract ab- 
sorption signals from bright GRB afterglows. Here we examine the 
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Figure 15. Cumulative distribution of equivalent width of the 21 cm ab- 
sorption lines around z ~ 10. The dot-dashed, short dashed, solid, dotted, 
and long dashed curves are computed assuming fx = 0, /x = 0.05, 
fx = 0.1, /x = 0.2, and fx = 1.0, respectively. 



observability of 2 1 cm signals on both spectra of these background 
sources. 

As for the upcoming and planned low frequency interferom- 
eters, a spectral resolution of 1 kHz is achievable for LOFAF0 
and SKA0. From the theoretical spectrum computed above, we get 
about one absorption line in every 8.4 kHz at z ~ 10 on average 
when we use M m in = 1O 6 M0 (line overlapping is accounted), and 
the line density (per observed frequency interval) is lower for lower 
redshifts or higher M m i n . On the other hand, the line width ranges 
mostly from ~ 1 kHz to ~ 5 kHz for halos of different masses. 
So the instruments can marginally resolve these 21 cm lines. While 
resolving the detailed line profile is probably out of reach, the line 
counting is feasible as long as sufficiently bright radio sources can 
be found at high redshift. 

To produce mock spectra matching real observations, we de- 
generate the theoretical spectrum to a resolution of 1 kHz, add 
Gaussian noise onto each pixel with the signal-to-noise ratio 
S/N — 5, and convolve it with a continuum of GRB afterglow 
as well as a radio spectrum of quasar. Two synthetic spectra are 
illustrated in Fig[T6] with the left panel for a GRB afterglow and 
the right panel for a quasar. We take fx = 0.1 for both cases. The 
flux density of GRB afterglow is scaled to 100/iJy at 200 MHz and 
z — 6, which is achievable for the afterglow of an energetic GRB 
produced by explosion of a massi ve metal-free star at high redshift, 
with isotropic energy of 10 erg dloka & Meszarosll2005h . and the 
spectral index is taken to be 2 (F v oc v 2 ), which is appropriate 
for the s ynchrotron self-absorption spectrum at the frequencies of 
interest jFrailll2003l) . The quasar flux density is scaled to 20mJy 
at 120 MHz and z — 10, with a spectral index of —1.05 as fit- 
ted to the radio sp ectrum of the powerful radio galaxy Cygnus A 
dCarilli et alj2002h . 

Apart from spectral resolution considerations, we also need 
high enough sensitivity to observe the absorption lines. In other 
words, with the planned instruments, the background source has to 
be bright enough to get the decrement of flux density higher than 
the detection limit. The minimum detectable flux density of an in- 
terferometer is related to the system temperature T sys , the effective 
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aperture area A e n, channel width Ai/ cn , integration time t[ n t, and 
the signal-to-noise ratio S/N by 



AF m 



2 k B r s . 



A eBy /Ap ch t in t N 



(33) 



The real signals of minihalos or dwarf galaxies are their addi- 
tional absorptions with respect to the absorption by the IGM, i.e. 
AF = F v exp(— tigm) — F v exp(— r). Equating this flux decre- 
ment to the detection limit, we get the minimum flux density of the 
background source required to observe the absorption lines: 



F min = 542 M Jy^^J 
/5000 ir^K" 1 



0.1 



V Acf[/T s}l 
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tint J 
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(34) 



where the ratio A c ff/T sys is an intrinsic parameter describing the 
sensitivity of an interferometry array, and we use the value for SKA 
here. For the GRB afterglow, the integration time is limited by its 
fading time scale. Typically, after a bright, short-lived radio "flare" 
at early times, the subsequent evolution of the radio afterglow can 
be described by a slow rise to maximum, followed by severa l seg- 
ment s of power-law decays with a timescale of ~ 100 days dFraill 
l2003h . Here we have assumed a reasonable integration time of 30 
days, and find that a minimum flux density of ~ 500 /iJy is re- 
quired to detect the absorption lines with the resolution of 1 kHz 
for the case of fx = 0.1. 

As seen from the continua (dashed lines) in Fig[l6] the flux 
density of our prototype GRB afterglow is more than one order of 
magnitude lower than this limit. Note that this is already an ener- 
getic GRB which is one order of magnitude brighter than normal 
GRBs. It may be possible that there are even brighter GRBs, but 
for most GRBs it seems that the radio afterglows would be too dim 
for being used as background source in observations with such a 
high spectral resolution. If we by chance find a quasar at very high 
redshift during the early stages of reionization, the signals could be 
easily detected. Especially, if one stacks together several lines to 
get an average profile, it will hopefully reveal the horn-like profiles 
we found. 

Alternatively, we may try broadband observations with lower 
resolution. In this case, we could use the standard measurement of 
Da, i.e. the mean (relative) flux decreme nt in each ban d, in analogy 
to the Lya forest experiments of quasars l lRauchll998h . We redefine 
the mean flux decrement with respect to the continuum after the 
absorption by the IGM: 



D A = ( /lGM ~ /obs ) = (1 

/iGM 



(35) 

where the angular brackets represent the average over each band 
and T e fl is the effective optical depth of the band. Da measures the 
excessive absorption by minihalos/dwarf galaxies compared to the 
diffuse IGM, and the different values of Da measured in differ- 
ent bands represent the absorptions at different redshifts, showing 
the evolution of non-linear structures during the epoch of reion- 
ization. In terms of Da, the flux decrement can be written as 
AF = F y exp(-riGM) - F v exp(-Tefi) = F v Da exp(-TiGM), 
then the requirement of the background source is 
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Figure 16. Synthetic spectra of 21 cm absorptions against a GRB afterglow (left panel) and a quasar (right panel) with 1 kHz resolution. The corresponding 
redshift is z ~ 10, and the unabsorbed continua are shown as the dashed (blue) lines, /x = 0.1 is assumed for both spectra. 



With the broadband observation, the major concern is that the 
available length of the line of sight is limited. On one hand, we 
wish to put the background source as far as possible in order to 
get enough bands to see varying Da values at different redshifts. 
On the other hand, as the source moves farther, it gets dimmer, 
lower the average signal to noise ratio, then an even broader band 
is required to detect the signal, thus beyond a certain distance the 
number of useful bands might decrease again. An optimal source 
redshift zgrb can be found to maximize the number of available 
bands. Considering the applicability of our model, we set the lower 
limit of the redshift to be zn m = 8, and find that the optimal red- 
shift for the GRB is zgrb ~ 9.8. This redshift depends mainly on 
the lower limit of the redshift we set, but of course it is also depen- 
dent somewhat on the reionization history, e.g. Da which evolves 
slowly with redshift. This does not depend on the observational pa- 
rameters. The optimal redshift is quite encouraging given that we 
have already seen 2 GRBs beyond z = 6. 

The existence of an X-ray background and its intensity are 
crucial in determining the observability of the signals. If there is 
no X-ray background, then with a GRB at redshift 9.8 and with the 
pixel resolution of Au c h = 1-38 MHz, we can get 19 bands, then 
19 values of Da- We plot the Da as a function of observed fre- 
quency in the upper panel of Fig[l7] and with the Gaussian noises 
of S/N = 5, the expected spectrum after the absorption by miniha- 
los and dwarf galaxies F v is shown in the bottom panel of Fig|17| 
The original continuum F c and the flux density absorbed by the ho- 
mogeneous IGM -Figm are also shown for comparison. With SKA, 
we could hopefully detect the non-zero Das if /x — or it is 
extremely small, which are signals from both minihalos and dwarf 
galaxies. In addition, Da decreases statistically with the decreas- 
ing redshift, which is a clear signature of the evolution of non- 
linear structures during the epoch of reionization. In the case of 
fx = 0.05, we could marginally observe 2 pixels on a spectrum 
from 131.5 MHz to 158 MHz, one of which has a band width of 
~ 8 MHz (131.5 - 139.5 MHz) with the mean flux decrement of 
D A ~ 0.017, and the other has a band width of ~ 18.5 MHz 
(139.5 - 158 MHz) with the mean flux decrement of Da « 0.009. 
However, if fx takes our fiducial value 0.1, then even one band 
spanning the whole spectrum from the source to the lower redshift 
limit is still not broad enough to make the observation feasible. 
The situation gets worse for higher values of fx. Thus, the 21 cm 
broadband observation against high redshift GRBs can be a power- 
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Figure 17. The mean flux decrement with respect to the continuum ab- 
sorbed by the diffuse IGM (upper panel) and the synthetic broad band spec- 
trum with 1.38 MHz resolution (bottom panel) in the case of fx = 0. The 
19 pixels show the frequency range corresponding to the redshift range from 
9.8 to 8. In the bottom panel, the dashed, dot-dashed, and solid lines are the 
original continuum flux density, the flux density absorbed by the homoge- 
neous IGM, and the expected flux density after the absorption by minihalos 
and dwarf galaxies, respectively. 



ful probe of the presence and the intensity of the early X-ray back- 
ground and the thermal evolution during the epoch of reionization. 



6 DISCUSSION 

We have modeled in detail the gas density and velocity profiles, 
ionization state, temperature profiles, and Lya photon production 
for both minihalos and dwarf galaxies during the early stages of 
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reionization when the IGM was still patchy. We also take into ac- 
count an early X-ray background, which could partially ionize and 
heat the IGM, and suppress the formation of low mass minihalos. 
Using the detailed model, we investigate the spin temperature of 
neutral hydrogen at different radii and the optical depth profiles of 
21 cm absorption lines, for various impact parameters, halo masses, 
and redshifts. 

We find that the Lya background and the Lya photons from 
recombinations in the slightly ionized IGM couple the spin temper- 
ature to the kinetic temperature outside minihalos or HII regions of 
dwarf galaxies, and the coupling is already strong at redshift 10. 
The Lya background photons are blocked at surfaces of halos, and 
the Lya photons from recombinations are negligible inside low- 
mass halos, though it dominates over the collisional coupling in 
high-mass halos in which the gas is partially ionized by collisions. 
The collisions also couple the spin temperature to the gas kinetic 
temperature effectively at the center, but this coupling decreases 
and the spin temperature decreases with increasing radius. 

The infalling velocity of gas around the minihalos/dwarf 
galaxies plays a very important role in determining the profiles of 
the 21 cm lines, double-peaked horn-like profiles are produced for 
a vast range of parameters. The line profile of a dwarf galaxy also 
depends on the size of its HII region, and hence its radiative proper- 
ties. The horn-like profile disappears for galaxies which are large, 
or with top-heavy IMF so that the HII region expands beyond the 
turning point of gas velocity. The optical depth of a dwarf galaxy 
is lower for lines of sight penetrating through its HII region, while 
a sufficiently large HII region will lead to an optical depth trough 
rather than an absorption line. 

With the line number density based on halo mass function 
and a physically motivated criterion for star formation, we gen- 
erate synthetic spectra of 21 cm forest by Monte Carlo procedure, 
and calculate the cumulative distribution of equivalent width of the 
absorption lines. Most of these lines have equivalent widths around 
0.03 ~ 0.3 kHz for an X-ray background intensity parameter of 
fx = 0.1, and the number of strong signals with large equivalent 
widths decreases significantly with increasing fx- We then study 
the observability of these signals. For high resolution (1 KHz) ob- 
servations, the GRB radio afterglows are too dim to be used as the 
background, but absorption lines should be easily detected for a 
high redshift quasar. It is exciting to know that the Pan- STARRS 
(the Panoramic Survey Telescope And Rapid Response System) is 
being developed which will be able to detect quasars up to redshift 
~ 7 and aims to find ~ 20 — 50 quasars at z ~ 7. For broadband 
observations, it is also possible to detect the absorptions against 
GRB radio afterglows if there is no X-ray background. Setting a 
lower redshift cut-off at zum = 8, we find that the optimal redshift 
for the GRB is zgrb ~ 9.8. With a sensitivity of SKA, a signal- 
to-noise ratio of 5, and a reasonable integration time, we could get 
measurements of mean flux decrement Da for 19 bands along the 
line of sight, each with a channel width of 1.38 MHz. In this way 
we could detect not only the signals from both minihalos and dwarf 
galaxies, but also their evolution during the epoch of reionization. 
However, the detectability of 21 cm signals is very sensitive to the 
presence of an early X-ray background. If an early X-ray back- 
ground existed but was not strong, taking fx = 0.05, we could 
marginally observe 2 pixels along the line of sight towards a GRB 
which is located at the optimal redshift of 9.8. Nonetheless, for the 
value of /x = 0.1 or higher, the signal will be impossible to be de- 
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tected. Therefore, the 21 cm absorption could be a powerful probe 
of the presence/intensity of the X-ray background and the thermal 
history in the early universe. However, we note that it is difficult to 
find a radio source before the IGM has been significantly ionized 
or heated, especially for a very bright high redshift quasar. 

IFurlanetto & Loebl J2002h studied the 21 cm forest signals of 
minihalos and early galaxies. Here, we have re-investigated this 
problem with different and more detailed modeling of various prop- 
erties of these nonlinear objects, an early X-ray background, as 
well as the Lya background during the epoch of reionization. We 
found stronger absorption signals from both minihalos and dwarf 
galaxies for an early X-ray background not higher than the level 
today. Given the many different model parameters adopted, this is 
not unexpected. There are many differences in the details of mod- 
eling between the two p apers, but the main diffe rence seems to 
be the IGM temperature. IFurlanetto & Loebl d2002l) have assumed 
a heated IGM with a simple form of its evolution with redshift, 
which is already heated up to 1000 K at z = 10. This lies be- 
tween our cases of fx = 1 and fx = 5. But in our fiducial 
model, the gas temperature in the IGM is about 35 K at the same 
redshift for fx = 0.1. This is appropriate for the early stages of 
stellar reionization , when the percolation has not occurred yet. In 
IFurlanetto & Loebl J2002h . the gas structure in the infalling region 
around minihalos/dwarf galaxies is modeled with a self- similar so- 
lution of secondary infall found by iBertschingerl 1 19851). w hile we 
have used the gas infall model developed bv lBarkanal d2004h which 
is based on the extended Press-Schechter model and spherical col- 
lapse. The Bertschinger's solution has a power law density profile 
of p oc r~ 2 ' 25 , which is much steeper than the Barkana's predic- 
tion. Our density and peculiar velocity structure of the infalling 
gas, with the lower spin temperature outside the minihalos/HII re- 
gions, produce a higher optical depth and a horn-like profile, which 
was not found in[F urlanetto & Loebl d2002h . For the dwarf galaxies, 
IFurlanetto & Loebl J2002h considered protogalactic disks, whereas 
we assume here a spherical symmetry in the gas density distribu- 
tion since the earliest galaxies are not likely to have large angular 
momentum. 

In addition to the early X-ray background, an obvious uncer- 
tainty in our model is the production of soft X-rays by dwarf galax- 
ies in the early universe. Especially with the normal IMF model-A 
and relatively high metallicity, the amount of soft X-rays emitted 
after a star burst is very uncertain dSchaererll2003l) . However, this 
amount of soft X-rays produced by stellar sources is always neg- 
ligible as compared to the background X-rays even if fx — 0.05. 
Another uncertainty comes from the gas density in minihalos and 
dwarf galaxies. We have neglected the possible change of density 
profile after the star formation. The gas density profile could be 
modified by the expansion of the HII region when the ionizing front 
changes nature from R-type to D-type. There are also uncertain- 
ties in the star formation history in the dwarf galaxy. We have as- 
sum ed that the feedback ef fects quench subsequent star formation 
(c.f. lOmukai&Nishilll999l) . so that a single star burst is produced. 
We have also assumed a smooth gas density distribution, neglect- 
ing cool dense gas clumps from which the first stars are likely to 
form. Thus we may have overestimated Run, and the optical depth 
could be slightly higher if they are accounted for. However, the gas 
fraction in minihalos or dwarf galaxies could also be lower than the 
cosmic mean value dNaoz et al.|[2009h . then the optical depth will 
be lower. Also, the results depend on the assumption of / csc , which 
is taken to be 0.07, in agreement with current observations. 

As illustrated in Fig|9]and Figll II minihalos and dwarf galax- 
ies exhibit distinct optical depth profiles mainly due to the different 
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ionization state and coupling physics, so they are potentially dis- 
tinguishable. Although we may not be able to resolve the line pro- 
files with the upcoming and planned instruments, it is encouraging 
to be able to distinguish their features in a statistical way, given 
that dwarf galaxies span a different halo population from miniha- 
los which cannot host stars. We reserve this investigation to future 
works. 
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